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Abstract 

We discuss several steady-state rotation and oscillation modes of the planar parametric rotator and 
pendulum with damping. We consider a general elliptic trajectory of the suspension point for both rotator 
and pendulum, for the latter at an arbitrary angle with gravity, with linear and circular trajectories 
as particular cases. We treat the damped, non-linear equation of motion of the parametric rotator 
and pendulum perturbatively for small parametric excitation and damping, although our perturbative 
approach can be extended to other regimes as well. Our treatment involves only ordinary second-order 
differential equations with constant coefficients, and provides numerically accurate perturbative solutions 
in terms of elementary functions. Some of the steady-state rotation and oscillation modes studied here 
have not been discussed in the previous literature. Other well-known ones, such as parametric resonance 
and the inverted pendulum, are extended to elliptic parametric excitation tilted with respect to gravity. 
The results presented here should be accessible to advanced undergraduates, and of interest to graduate 
students and specialists in the field of non-linear mechanics. 

1 Introduction 

A physical system is said to be parametrically excited when some of its parameters vary periodically with 
time. One of the simplest examples of a parametrically excited system consists of a rigid rotator (or, in the 
presence of gravity, a physical pendulum) whose suspension point follows a periodic trajectory, with both 
rotator and suspension point restricted to move on a plane. The basic modes of motion of that system are 
rotation and oscillation. Depending on the type of parametric excitation, its amplitude, and the strength 
of damping, many different types of rotation and oscillation modes may be encountered in its spectrum of 
steady states. Furthermore, those basic modes combine to form complex motions in which rotation and 
oscillation alternate in regular or irregular patterns, including chaotic motion. A qualitative description of a 
variety of different modes of motion in a parametric pendulum, and some further references, is given in the 
introduction of [T]. 

The parametric rotator does not seem to have received much attention in the literature, despite the 
fact that its rich non-linear dynamics (see, for example, [2]) is interesting for both research and education 
purposes. The parametric pendulum is discussed in the literature in relation to its two most remarkable 
dynamical phenomena. One of them is parametric resonance, that occurs in systems possessing normal modes 
of oscillation when the period of the parametric excitation is an integer or half-integer multiple of the period 
of a normal mode. Parametric resonance in a pendulum with vertically or horizontally moving suspension 
point is discussed in the textbooks 0|4], and in [SHH]. Parametric resonance has also been discussed in the 
context of the elastic or spring pendulum [SHU] , a linear oscillator jlS^JJj , elastic strings (TUllin] , electronic 
circuits [TUlfTB] , and many other systems. 

The other striking dynamical aspect of the parametric pendulum is the so-called "inverted pendulum." 
A vertical parametric excitation, with sufficiently large amplitude and frequency, causes the upright, usually 
unstable, equilibrium position of the pendulum to become stable. The dynamical stabilisation involved in 
this phenomenon is of interest by itself, and also because it constitutes a mechanical analogy for similar 
mechanisms in more complex systems |17[|18| . A detailed study of the inverted pendulum was given by 
Kapitza [TH] . Textbook treatments are given in [2U[[^ . In the more recent literature a significant number of 
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studies of the inverted pendulum exist, both experimental P"51I^2H^ and theoretical [T]|^]|51 IT71I^T]|25H5^ (see 
also the references cited in those papers) . Most of those theoretical works deal with the non- linear equation of 
motion of the parametric pendulum by linearising it in the small-oscillation approximation |21ll25ll26ll29ll31j in 
which it reduces to a Mathieu equation (see, e.g., |33j and references therein), or within the framework of the 
effective potential by separation of slow and fast variables ,1,2^,20,, or by considering non-sinusoidal motion 
of the suspension point [5l|27ll28]. In most cases damping is justifiably neglected. All of those theoretical 
approaches lead to important insights into the dynamics of the inverted pendulum, and to quantitative 
results such as, for example, the determination of its region of stability. 

In this paper we adopt a complementary point of view. We look for steady-state solutions to the 
damped, non-linear equation of motion of the parametric rotator and pendulum. Our approach is based 
on the observation that the basic steady-state rotation and oscillation modes are mathematically simple, 
in the sense that they can be expressed in terms of elementary functions in some perturbative expansion. 
For such basic modes the non-linearity of the equation of motion manifests itself in the multiplicity of 
different regimes, characterised by a few dimensionless parameters, leading to different types of solution 
and requiring different perturbative expansions. Here, we restrict ourselves to small parametric excitations 
and low to moderate damping (both to be quantified below), and set up a perturbative expansion in the 
small excitation amplitude. In that regime the spectrum of basic modes of motion is simpler than at larger 
amplitudes, though still complex, but more interesting than in the limit of high friction. The cases of large 
parametric excitations and/or large friction could, in principle, be handled by straightforward modifications 
of our perturbative method. Some of the steady-state rotation and oscillation modes studied here have not 
been discussed in the previous literature. Other well-known ones, such as parametric resonance and the 
inverted pendulum, are extended to elliptic parametric excitation tilted with respect to gravity. 

A pedagogical advantage of the approach adopted here is that it yields numerically accurate approximate 
solutions to the equations of motion of realistic non-linear systems like the parametric rotator and pendulum 
by elementary mathematical methods. The only technical requirements are the ability to perform power- 
series expansions and to solve inhomogeneous ordinary second-order differential equations with constant 
coefficients, which should be familiar from a course on analytical mechanics. Thus, we hope that the 
analysis in the following sections may contribute to our understanding of non-linear mechanics, and to make 
the latter more widely accessible to both students and educators. 

The paper is organised as follows. In the next section we derive the equation of motion for the rotator 
with general elliptic motion of its suspension point, subject to gravity and damping, and establish our 
notation and conventions. In section [3] we discuss the two-dimensional rotator with small linear parametric 
excitation, in the absence of gravity, explaining in detail the perturbative method to solve the equation of 
motion for the steady-state of the system. The same perturbative approach is applied in the subsequent 
sections. The extension of the results of section [3] to general elliptic parametric excitation is carried out in 
section m The inclusion of gravity is considered in section [SJ where we study rotation and oscillation modes 
of the parametric pendulum elliptically excited at an arbitrary angle with gravity. Steady states related to 
parametric resonance are discussed separately in section [6l Finally, in section [7] we give our summary and 
final remarks. 

2 The parametric rotator 

We consider a two-dimensional rotator, that is, a rigid body free to rotate about one of its points whose 
trajectory in space is fixed. We call the rotator symmetric if its suspension point is at its centre of mass, 
otherwise asymmetric. If a rotator is symmetric, a spatially constant force field, such as an inertial or 
gravitational force, cannot exert a torque on it. To be specific, we will think of a thin homogeneous bar of 
mass M and length L with its suspension point at a distance Li from one of its ends and L2 from the other 
(I/i + L2 — L, Li > ^2)- If we denote by tq the position of the suspension point, and by ^ a coordinate 
along the bar, —L2 < £, < Li, the points on the bar are parameterised by r(^, t) = ^ sin 0a; — ^cos^z -1- tq. 
It will be convenient to choose a coordinate system such that the gravitational force points not in the 
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direction — z, but forming an angle —a with it. The rotator's Lagrangian is then, up to total derivatives, 

1 /"^^ M ■ f^^ M 

L = — / — r(^,t)^— / — (7(cosq;zo — cosacos^^ + sinaxo + sinasin^^) 

2 J-L2 L J-L2 L 

= ^le^ + ^M{Li - L2) (sin 61x0 - cos 6*^0) + ^Mfi^ + ^Mg{Li - L2) cos{9 + a) 
— Mg{cos a zq + sin axo) , 

where T = M/Z{L\ + L2 — L1L2) is the moment of inertia with respect to the suspension point, and where 
in the second equality we dropped a total derivative. In what follows we always take into account the effect 
on the rotator of friction, in the form of viscous damping proportional to its angular velocity. Adding such 
term to the equation of motion for from the Lagrangian ([T]) we obtain 

l(e + \e) = ^M{Li - L2) (cos^^o + sin6lzo) - ^Mg{Li - L2)sm{e + a) , (2) 

A being the damping coefficient. On the left-hand side of ([2]) the first term corresponds to the rotator's 
inertia and the second one to damping due to friction. On the right-hand side, the first term describes the 
torque applied on the rotator by the inertial forces due to the accelerated motion of its suspension point, and 
the second one the torque due to gravity. Notice that the right-hand side vanishes for a symmetric rotator. 

We restrict ourselves to the case where the suspension point traces an elliptic trajectory about the origin 
with a single frequency 7, 

Xo{t) = -ro sinesin(7t) , zo{t) = tq cos(7t) , (3) 

with To the length of its semi-major axis, tq sine that of its semi-minor axis. Notice that we have chosen 
the z axis along the major axis of the ellipse ([3]), which explains why we took gravity along an arbitrary 
direction in ([l} . Substituting ([3]) in ([2]) we obtain the equation of motion for the rotator with an oscillating 
suspension point. It is convenient, however, to change the time variable in ^ to dimensionless time r = jt, 
so as to express the equation in terms of a minimal set of dimensionless parameters. After doing so, we 
obtain, 

^ + + A cos -f A f^lil) sin ^ + r sin(0 + a) = , (4) 
ctr^ dr To ro 

with A = A/7 the dimensionless damping coefficient and 

_ MjLl ~ Ll) ro M(Ll-Ll) g 

21 L ' 2Z L72 ■ ^ ' 

Equation ^ is generally valid for any rotator with its suspension point oscillating as in , and coefficients 
A, r depending on its geometry and mass distribution. For a symmetric rotator A = = F. For the 
homogeneous thin bar discussed above, A and F are given by ([5])- 



3 Linearly excited parametric rotator 

We consider first the case in which the plane of the rotator is horizontal, so that gravity does not play a role 
(F = in ([4])), and the motion of its suspension point is linear (e = in (|3])). The equation of motion Q 
then takes the form 

Notice that in ([B]) we can always choose the time origin so that A > 0, and that the equation is invariant 
under 6 —9. We restrict ourselves to the regime of small parametric excitations < A < 1. 

3.1 Basic solution 

We look for steady-state solutions to ([6|) with 9(t) linearly rising with t, and possibly performing small 
oscillations about that linear trajectory. Thus, we are led to propose a solution of the form 

0(r) =c^T + eo + /(r)A , w^O, (7) 
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with Qq a constant phase. It is physicahy reasonable to include such constant term in 9(t) because, due 
to the motion of the suspension point, the z axis is a privileged direction in space, making ^ sensitive to 
the phase of 9. For instance, 9 — m: with integer n are equilibrium solutions to ©, but any other constant 
is not. The last term in ([7]) represents small oscillations with an amplitude of 0(A), which requires the 
unknown function /(r) to be bounded for < t < cxd. 

If we substitute (O in ([5]) and let A — 0, we find that luX — > 0. We must therefore assume that A — 0(A), 
which expresses the physically reasonable condition that, for small parametric excitations to cause steady- 
state rotations, friction must be correspondingly small. (The relation between A and A is made precise 
below in equation (fTTj).') Thus, substituting ([7]) in ([6]) and neglecting terms of 0(A^) yields 



= —uj^ — i sin {{uj + 1)t) cos 0o ^ ^ cos ((cj + 1)t) sin 6o 
dr-' A 2 2 

— — sin {{oj — l)r) cos Qq — — cos ((cj — 1)t) sin Qq 



(8) 



where A/A remains finite as A — > 0. Equation ([S]) is inconsistent with our assumptions, since the constant 
term on its right-hand side leads to an unbounded perturbation /(r) = — 1/2(wA/A)t^ + growing 
quadratically with r. The only way around this inconsistency is to set lo = ±1, to obtain, 

= T^-^sin(eo)T^sin(2T±eo) , = ±1 . (9) 

Setting now sin(0o) = T2A/A, the r-independent terms on the right-hand side of ([9]) vanish, and the 
equation admits a bounded solution /(r). We have then a solution to ([6]) of the form 9{t) — ±9o{t) + 0(A) 
with, ^ 

6'o(r)=r + eo, 60 = - arcsin(2A/A) + 2n7r . (10) 

This result shows that for small A there cannot be parametrically excited steady-state rotations unless 

A<A/2. (11) 

We notice also that the steady-state phase 0o is completely determined by the non-linear dynamics and 
the damping, without reference to the initial conditions which, for times t 3> 1/A, are completely erased 
by the damping interactions. We conclude that, in physical units of time, at small A and A satisfying 
(jlip . steady-state rotations occur only with the same angular frequency 7 as the parametric excitation, 
9{t) = ±(7< + Oq) + C(A), either clockwise or counter-clockwise. 



3.2 Perturbation theory: steady-state rotation 

We now develop the considerations of the previous section into a systematic perturbative approach. We 
assume a solution to ^ of the form, 

°° A" 

9{T) = 9oir) + Y,^ifn{r) + en) , (12) 
n! 

n— 1 

with 9q{t) given by pUj) . We consider only positive-frequency solutions 9{t) since, as noticed above, negative- 
frequency ones are just given by —9{t). The perturbations fn{T) in (|12p are assumed to be bounded in 
< r < 00, uniformly with respect to n. We have explicitly separated from /„ a constant term 0„, so we 
further require that fn{T) must not contain constant terms. The constants 0„ are perturbative corrections 
to the zeroth-order phase 0o in (fTU]) . 

Replacing in ^ and expanding to first order in A leads to equation ^ (with upper signs). With 
Oo given by (jlOp . and the condition that /i(t) must not contain constant terms, we get a unique solution 
to (HI) 

/i(t) = isin(2T + eo) . (13) 

In order to have a complete first-order solution we still need to find the constant 0i, which is fixed by 
requiring consistency of the second-order equation. 
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To second order in A equation ^ yields, 

Ia^^ + ~XA^ + Ia^ cos(2r + Qo){fi{r) + Q,) + ^A^ cos(eo)(/i(r) + Oi) = . 
Substituting fi in this equation we get, after some rearrangements, 

^ = cos(eo)ei - 4 '^in(2r) + ^ sin(2r + 260) - sin(4T + 260) • (14) 
ar^ 16 16 16 

The first term on the right-hand side is a constant that must vanish for /2 to be bounded, therefore 

61=0. (15) 

We can now find /2 by twice integrating (jl4p with respect to r. The requirement that /2 should not contain 
zero-frequency terms leads to a unique solution, 

/2(t) = ^ sin(4r + 260) + ^ sin(2T) - 1 sin(2r + 260) . (16) 

With this result the solution to ([6|) to 0(A^) is determined up to the constant 82, which results from 
requiring consistency of the third-order equation. 

We proceed further to consider the terms of 0(A'^) in equation © with (IT^ . After substituting /i and 
/2 in the third-order equation, somewhat lengthy algebra yields 

d/^ f / 3 15 \ 3 27 87 

" ( - 2 cos(0o)e2 + ^ sin(eo) j - - cos(2r + 60)82 - ^ sin(2r - Bo) + sin(2r + Go) 

3 45 15 9 ^^"^^ 

- ^ sin(2r + 360) - ^ sin(4r + 60) + ^ sin(4r + 360) - sin(6T + 360) • 

As before, boundedness of means that the constant term in brackets must vanish, so that, 

e2-^tan(eo). (18) 

Now equation (jl7p can be integrated, taking into account that /3(t) must not contain constant terms, to 
yield, 

3 27 87 3 

/3(t) = -62 cos(2t + Go) + ^ sin(2T - 60) - ^ sin(2T + Q„) + — sin(2T -f 360) 

45 15 1 

+ ^ sin(4T + eo) - ^ sin(4r + 360) + sin(6r + SOo) 

The third-order correction is completed with the equality 



(19) 



e3 = i^sin(eo) , (20) 

required for consistency of the fourth-order equation. 

The above results constitute a perturbative solution to ([5]) as a series expansion in A, valid in the regime 
< A < 1 and when the condition (|lip holds. The zeroth-order solution is (fTO|) . The first, second and third 
order solutions are given by, 

0i(t) = 0o(t) + /i(t)A , 02(T) = 0l(r) + (e2 + /2(T))^ , ^3 (t) = ^2 (t) + (63 + /s (t)) ^ . (21) 

Several remarks about this perturbative solution are in order, (a) In the perturbative expansion ()12p we 
could have included corrections to the angular frequency, with a first term in ()12|) of the form (1 -|- wiA -|- 
a;2 A^/2 + ■ ■ ■)t instead of just r. All the corrections a;„, however, turn out to vanish, (b) The results 
(fT6| . (|19p. for /i,2.3 suggest that the perturbative expansion parameter is actually A/8 rather than A, so 
the domain of applicability of perturbation theory should be larger than expected. This is confirmed below 
in section 13.51 on numerical results, (c) Qq is determined only up to 2n7r. The additional n turns typically 
build up during the transient state. Our perturbative solutions hold valid only for times t ^ 1/A after all 
transients have died off, and depend on A only through the phase Qq. (d) The perturbations fn{T) contain 
only even frequencies (in physical units, only even multiples of 7). Thus, 9{t) — t is periodic with period tt. 

These results are numerically verified below in section 13. 5[ where they are compared with numerical 
solutions of ([6]) and found to be in tight agreement with them, each order significantly improving on the 
previous one. Those verifications justify our perturbative procedure, together with the implicit hypothesis 
that the dependence of 9{t) on A is analytic at A = 0. 



5 



3.3 Steady-state oscillations 



When condition (jlll) is not satisfied the steady-state rotation described in the previous section does not exist. 
In that case, the rotator may reach the equihbrium solutions 9 = mr. There is another type of steady-state 
solution that is possible in that case, however, describing oscillations about (2n + l)7r/2. To obtain those 
solutions we assume they have the form, 

°° A" 

=0o + 5]^(/n(r) + e„) . (22) 

n=l 

Substituting this form into the equation of motion ^ leads to perturbative equations for /„ that can 
be solved by the same method as in section 13.21 We will omit the details of the procedure for brevity, 
commenting only on the form of the solution. 

The condition that the perturbation /2(t) be bounded leads to sin(29o) — 0, whereas the equations at 
third, fourth and fifth order yield Qi — Q2 — = 0. If we set Qq — /ctt, for some integer fc, we recover 
the equilibrium solutions and all corrections 0„, /„ must vanish because Qq is itself an exact solution to 
the equation of motion. Another possible choice is Qq — {2k + l)7r/2, with k an integer. In that case the 
solution is non trivial. At first order we obtain, 

6*1 (r) = 60 + A sin(eo) cos(t) . (23) 

The equation for /2 yields the second-order solution, 

6»2 (r) = 9i (r) - AA sin(eo) sin(r) , (24) 

and at third order we find 

Osir) = 02{t) - sin(eo) (^^A^ cos(r) + Ap cos(t) + ^A^ cos(3t) j . (25) 

The perturbative method does not yield information on the domain of validity of these solutions in parameter 
space. Numerical analysis of ([6]) suggests that, in the perturbative regime A, A < 1, the steady state described 
by (123])-(I25]) exists only for A < A. 



3.4 Order-of-magnitude estimates 

The quantitative meaning of the parameter A in ([6]) is apparent from its definition ([5]) for a homogeneous thin 
bar, or its correspondingly modified expressions for other rigid bodies. Thus, for example, for a homogeneous 
thin bar with Li = 5.4 cm, L2 = 4.6 cm and excitation amplitude ro = 0.5 cm we have A = 0.016. For 
a larger, more asymmetric thin bar with Li — 13 cm, L2 = 7 cm and amplitude vq ~ 2 cm we obtain 
A — 0.094 ~ 0.1. By reducing the asymmetry or the amplitude we can get values of A smaller by several 
orders of magnitude. With larger asymmetries and amplitudes A can be made larger, A > 1. For values 
A > 2 the parametric-excitation amplitude must be larger than the length of the rotator itself, so rather 
than parametrically excited, the rotator should in that case more properjy be referred to as orbitally excited. 

In order to obtain some quantitative understanding of the parameter A in © , we consider the steady-state 
rotations described in section [3^ Once the steady state has been reached, we may cease the parametric 
excitation and find out how many turns the rotator completes before it stops moving. Since for a steady-state 
rotation of the form ([T2l) the (dimensionless) angular velocity is 1 -|- 0{A), the number of turns sought for 
must be approximately l/(27rA) independently of A. Indeed, after the motion of the suspension point has 
stopped, we find by numerically solving ([6]) with A = that the rotator performs slightly more than 15 
turns if A = 10"^, about 150 turns if A = 10"^, and about 1500 turns for A = 10"'*. From our everyday 
experience we know that the first case corresponds to the level of damping usually encountered in common 
household appliances such as fans or mixers. The second case corresponds to a much lower level of damping, 
found in more special mechanical devices such as some microprocessor fans. In all three cases we find that 
the perturbative results of section \3l2\ for the steady state agree with exact numerical solutions to ([6]) with 
high accuracy, as discussed in the following section. For much lower damping, say A < 10~^, which can be 
achieved in laboratory experiments, the steady state takes a long time to be reached. In those cases, it may 
be appropriate to modify our perturbative expansion by setting A = 0(A") with n > 1, or just solve (jB)) 
numerically for the entire time evolution, including the long transient state. 
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3.5 Numerical results 



We now turn to comparing numerical solutions to ([B]) with the perturbative ones found in the previous 
sections. We numerically solved ^ with standard commercial software |34j over an interval < r < Tmax 
with Tmax 3> 1/A, with appropriate initial conditions. The type of steady state reached by the system is 
insensitive to 6{0) in most cases, but may be sensitive to dO/dT{Q) depending on the values of A and A. 

In figure [T] we show steady-state rotations for three different values of A and A. With A = 0.1 and 
2A/A = 0.2, as shown in figure [1] (a) , the first-order approximation 0i{t) from (pij) already agrees almost 
exactly with the numerical solution. For these parameter values we numerically find steady-state rotations 
for any value of 0(0), and for 0.64 < d6/dT{0) < 2.13. Outside of that range of initial angular velocities the 
steady states we obtain numerically are the equilibria 9{t) = mr. 

In figure [1] (b) we set A — 0.8 and 2A/A = 0.025. The first-order approximation is very close to 
the numerical solution, and the second-order one 02 (t) is indistinguishable from it. We numerically find 
steady-state rotations for any 0{Q) ^ 2mT, and for all values of 10~^ < d9/dT{0) < 10^ we tried. 

Numerical analysis of ^ shows that the steady-state rotations described in section 13.21 for A < 1 exist 
also for larger values of A, although the perturbative expansion discussed there is not directly applicable 
when A ^> 1. (In the latter case, an expansion in powers of 1/A would be appropriate.) Nevertheless, 
as discussed at the end of section 13. 2[ the form of the expansion suggests that it may converge even for 
A > 1. This is shown in figured] (c), where A = 2.3 and 2A/A = 0.61. As seen in the figure, the agreement 
between the numerical and the perturbative solutions is very accurate at third order. Also, each successive 
perturbative order gives a significant improvement over the previous one, signalling good convergence of 
perturbation theory even at such large values of A (as long as damping is also sufficiently large). In this 
case we also find steady-state rotations for any 9{0) ^ 2mr, and for all values of d9/dT{0) over many orders 
of magnitude. 

The steady-state rotations described in section [3?2] and in figures [T] (a), (b), (c), exist for large open sets 
in parameter space, even in the perturbative regime A, A < 1, and are reached by the system for relatively 
large open sets of initial conditions, as described above. The steady-state oscillations described in section 
13.31 by contrast, are found numerically only for A <^ A. Under those conditions, the angular velocity quickly 
reaches its typical steady-state magnitude and, as a result, the mean value of 9 evolves very slowly towards 
its steady-state value and numerical solutions are computationally expensive. Also, unlike steady rotations, 
these oscillating steady-states are not robust in the sense that even small variations in the initial velocity 
may lead the system away from them and towards an equilibrium solution. An example of steady-state 
oscillations about 7r/2 is shown in figure [T] (d). In that case the numerical result is well described by the 
first-order solution (f23|) and, in the figure, indistinguishable from the second-order one ([24]) . 



4 Elliptically excited parametric rotator 

The simplest implementation of a parametric rotator is probably the well-known toy model described in |35j 
as a "rotator on a notched stick." In that case the trajectory of the suspension is well approximated by a 
highly eccentric ellipse. The elliptical excitation of the parametric rotator is of interest by itself, and also 
because it continuously connect the cases of linear (e = 0) and circular (e = ±7r/2) excitations, that lead to 
rather different types of steady states. In this section we consider the extension of the results of section [3] 
to the case of elliptic motion of the suspension point. 

The equation of motion (j4|), without gravity, takes the form, 

3-| + + ^(1 - sine) sin(6l + r) + ^(1 + sine) sin(6l - t) = , (26) 
dr-^ dT 2 2 

where —tt/2 < e < tt/2. For e = the suspension point oscillates along a line segment, and we recover ([6]). 
For < e < 7r/2 the suspension point moves counterclockwise along an elliptic orbit, which becomes circular 
at e = tt/2. Similarly, negative values ~tt/2 < e < describe clockwise motion of the suspension point. In 
this paper we restrict ourselves mostly to eccentric elliptic excitations, |e| < 1/2 that are close to the linear 
case considered in previous sections. Circular, or slightly-eccentric elliptic, motions of the suspension point 
(1 ^ kl < T^/"^) lead to qualitatively different solutions that will be discussed elsewhere |36) . 
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Figure 1: Numerical solution (black solid line) to equation ([6]), compared with perturbative solutions at 
zeroth (grey dashed line), first (grey dotted line), second (black dashed line) and third (black dotted line) 
orders as given in section |31 Steady-state rotations (see section [221 (a) ^ = 0.1, A = 0.01, 6(0) — 0, 
de/dT{Q) = 0.7; (b) A = 0.8, A = 0.01, 6i(0) = 0, de/driO) = 0.001; (c] A = 2.3, A = 0.7, 6i(0) = -0.1, 
d9/dT{0) = 1. Steady-state oscillations (see section [SS]) : (d) A = 0.001, A = 0.1, e'(O) = tt/2, d9/dT{0) = 0. 
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Equation (I^Bl) is invariant under the transformation — > —6*, e — ^ — e. Thus, if 6'(t, e) is a solution to 
for some e, then —0{t, — e) is also a solution for the same e. There are no equilibrium solutions to 
when e 7^ 0. When e = 0, the linear case, the equilibria are 0{t) — rnr, with integer n. 



4.1 Steady-state rotations 

We now look for counter-clockwise steady-state rotations of the form p2| . Let 0(t, eo) be such a solution 
for e = eo > 0. Since both the rotator and suspension point move in the same sense, we will say that 
9{t, Co) describes a "direct motion" of the rotator. By symmetry, —9{t, —cq) is a solution of ([26]) with 
e = eo, describing a clockwise motion of the rotator with a counter-clockwise moving suspension point, what 
we call a "contrarian motion" of the rotator. Similarly, 9{t, — eo) and —9{t, cq) are solutions to with 
e — —Co < 0, describing contrarian and direct motion, respectively. 

After substituting (fT2|) in p6| , a straightforward application of the perturbative method of section [3] 
leads to the perturbative solution we seek. At lowest order we get 

9o{T)=T + eo, eo = - arcsin ( — -A — ^ 1 + ^'^'^ ■ (27) 

\ 1 + sui e A / 

Therefore, in this case steady-state rotations occur only if 

A<^(l + sine), (28) 

which generalises condition (fTTj) to elliptic excitations of the rotator. We see also that contrarian motions 
are allowed by (j28p . for small enough A, except for e = — 7r/2. Thus, for circular parametric excitation there 
can be no contrarian motion. 
At first order we get, 

^i(r) = 0o(T)+A/i(r) , /i(t) = i(l-sine)sin(2r + eo) . (29) 

In particular, 0i = as in the linear-excitation case of section [3l The second-order solution takes the form, 

A\„ „ „ 5 

(1 — sine)' tan(Wo) , 

(30) 



^2(t) = 0i(t) + — (e2 + /2(T)) , 92 = ^(1- Sin e)2tan(eo) 



/2(t) - ^(cose)2sin(2T) - ^(cos e)^ sin(2r + 260) + ^(1 - sin e)^ sin(4r + 260) . 
Finally, at third order we obtain, 

^3(t) =02(T) + ^(e3 + /3(r)) , 63-^(1- Sin e)(cose)2sin(eo) , 

3 27 
/3(t) - -62(1 - sine) cos(2r + 60) + ^(cose)2(l + sine) sin(2T - 60) 

3 3 
+ ^(29(sine)^ -t- 9(sine)^ - 9sine - 29) sin(2T + Go) + ^ (cos e)^ (1-1- sine) sin(2r + SOo) (31) 

45 15 
+ ^(cose)^(l - sine) sin(4r-|-eo) - ^(cose)^(l - sine) sin(4r + SOq) 

+ ^ (1 - sin e)^ sin(6T + 360) . 

Some remarks on these expressions are in order. First, at e = the perturbative solution described above 
reduces to the one found in section |3^ for the case of linear excitation. Second, at e = 7r/2 the perturbative 
corrections 82,3 and /i,2,3(''') vanish. In fact, all 0„, /^(t), n > 1, must vanish in that case because for 
circular excitation 9o{t), as given by ()27|) . is an exact solution to (|26|) . Third, only even angular frequencies 
enter the perturbative corrections. 
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4.2 Steady-state oscillations 



We now look for solutions to (1^^ of the form (1^^ . Assuming A, A < 1, and substituting in (I^Bl) . leads 
us to a perturbative analysis completely analogous to that of section 13. 3[ so we omit the details for brevity. 
From (1221) we have the zeroth-order solution 



Oo{t) = eo, (32) 
with Oo undetermined at this order. The first-order equation yields, 

/(g)i(T) = -^(1 + sine)sin(T - Go) + ^(1 - sine)sin(r + Oq) , (33) 

Consistency of the second-order equation requires cosesin(26o) — 0, therefore, 

eo = ^fc7r, if cose 7^0, (34) 

for some integer k. In the case of circular excitation, cose — 0, Qq as well as all 0„ remain undetermined 
due to the invariance of (j26p under a one-parameter continuous symmetry. The details of that particular 
case will be discussed elsewhere |36j. With Qq given by (l34t . the equation for /2 leads to, 

/2(t) = i&i - A/A)(l + sine) cos(t - Go) + (61 + X/A){1 - sine) cos(t + Go) 

1 1 (35) 

- —(1 + sine)^ sin(2r - 260) + —(1 - sine)^ sin(2T + 260) • 
lb Id 

The constant 0i, undetermined at this order, is fixed by the consistency requirement for the third-order 
equation 

©i = X7^^c°«(2eo) . (36) 
A (cose)^ 

With this value for 81 we obtain a bounded solution /3(t) to the third-order equation. We omit the 
expression for /3(t), however, because it is too lengthy to transcribe here. Finally, 

62-0 (37) 

is required for consistency of the fourth-order equation. 

For e = these solutions must satisfy Thus, setting k — 2n in leads to 60 — nn, and we recover 
the equilibrium solutions of the linear case. In particular, all corrections 0^, fi, i > 1 must vanish, since 
Qq is an exact solution to (jH) in this case. It is easily seen that this indeed happens with the corrections 
(p3| and (|35 |) - (|37l) . On the other hand, by setting fc = 2n + 1 in (|34| we recover the oscillating solutions 
discussed in section [3731 



4.3 Numerical results 

The perturbative results for elliptic motion of the suspension point are compared to numerical solutions to 
I^E^ in figure[21 Figure[2](a) shows a contrarian steady-state rotation, with e = —1/2 and A = 0.1. Because 
the oscillation amplitude is small we plotted 9{t) — t instead of 0{t). The zeroth-order solution (j27l) is then 
shown in the figure as 6o{t) — t = Qq ~ —0.876. As also seen there, the first-order result (|29|) is very close 
to the numerical solution, and the second-order result (|30p agrees exactly with it. 

For the contrarian motion with 6(0) = shown in figure [2] (a) we find steady rotations for 0.987 < 
d9/dT{0) < 1.090. For other values of d9/dT{0) the steady state is oscillatory, about mr for some integer 
n. Although not shown in the figure, we mention here that for 9(0) — t:/8 and 7r/4 we could not find 
steady-state rotations, any d9/dT{0) leading to steady oscillations about mr instead. If we set A and A to 
the same values as in figure [2] (a), but e = -1-1/2, then with 9{0) = steady-state rotations are reached 
for 0.559 < d9/dT{0) < 1.519, with 61(0) = tt/B for 0.579 < d9/dT{0) < 1.483, and with 61(0) = 7r/4 for 
0.624 < d9/dT{0) < 1.419. These numerical results suggest that the interval of initial angular velocities 
leading to steady-state rotations decreases in length with increasing 9{0). Also, as intuitively expected, for 
fixed 9{0) that length is larger for direct than for contrarian motion. 
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Figure 2: Numerical solution (black solid line) to equation ([26| . compared with perturbative solutions at 
zeroth (grey dashed line), first (grey dotted line), second (black dashed line) and third (black dotted line) 
orders as given in section |4l Steady-state rotations (see section |4J|) : (a) e — —0.5, A ~ 0.1, A = 0.02, 
61(0) = 0,d6i/dr(0) = 1; (b) e = 7r/8, A = 1.2, A = 0.05, 61(0) = 0, d6l/dr(0) = 1. Steady-state oscUlations 
(see sectiongJl): (c) e = 7r/4, A = 0.1, A = 0.1, 6i(0) = 0, d9/dT{0) = 1; (d) e = 7r/6, A = 0.001, A = 0.1, 
e{0) = 7r/2, d9/dT{0) = -1.33 x 10"'*. 

In figure[2](b) we show a direct steady-state rotation, with e — tt/8 and A — 1.2. In this case the zeroth- 
order result ([27l) is given by Qq — —6 x 10~^. Like in the case of linear excitation, we see that perturbation 
theory converges for A > 1. As seen in the figure, the second-order result ([50)1 accurately reproduces the 
numerical solution, and the third-order result (j3ip cannot be distinguished from the numerical one. For these 
values of the parameters we found steady-state rotations for any value of 9{0) and for any < d9/dT{0) < 10^. 

We show a steady-state oscillating solution in figured] (c), with e = 7r/4, A = 0.1 — A. In this case 
condition (|28p is violated, so steady rotations are not possible. The value of Gq = 57r ~ 15.708 receives 
corrections A8i, A'^/693, so that the solution oscillates about 8 ~ 15.851. The accuracy of the perturbative 
results from section 14.21 is apparent in the figure, where no departure of the third-order result from the 
numerical solution is visible. 

In the case of linear excitation (e = 0), as discussed in section 1331 steady-state solutions oscillating about 
(2rt + l)7r/2 could be found (see figure [T] (d)), although finding them required some degree of fine-tuning of 
initial conditions. For e ^ those solutions seem to be unstable. An example is shown in figure [2] (d), with 
e = 7r/6, where the solution initially oscillates about Q = tt/2 but moves away from that value to end up 
oscillating about Q — n. For smaller values of e the growth of the mean value of 9{t) from 7r/2 is slower, 
but even for e as low as lO^'' we could not find numerically a stable solution of this type. 
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5 Elliptically excited parametric pendulum 



In the general case of elliptic motion of the suspension point, with the z axis chosen along the major axis of 
the ellipse and with the direction of gravity at an angle tt + a with that axis, equation Q takes the form 

(^Q ~ d6 A A 

— r + A— + — (l-sine)sin(e' + r) + — (l + sine)sin(6l-r)+rsin(6' + a) = . (38) 
dr^ dr 2 2 

This equation is invariant under the transformation 9 — > —9, e — e, a — > —a. Thus, if 9{t, e, a) is 
a solution to (1381) for some e, a, then —9{t, — e, —a) is also a solution for the same e, a. There are no 
equilibrium solutions to (j38p when either e 7^ or a 7^ 0. For e = = a the equilibria are 9{t) — mr, with 
integer n, as in the case F = 0. 

When A, A, F < 1 we can solve (ISS)) perturbatively for the steady-state solution along the same lines as 
in the previous sections. The parameter F defined in ([S]) may be made arbitrarily small by tilting the plane 
of the rotator with respect to the vertical, so that g itself is small; or by making 7 large, 7 > y/g/L; or 
by making small the rotator asymmetry ((Li — L2)/L in the case of the thin homogeneous bar discussed in 
section [2]). For arbitrary, but small, values of the parameters, we discuss below steady-state rotations and 
two different oscillation regimes. Furthermore, for F ~ 1/4, there are two additional types of solutions that 
we discuss in section El 



5.1 Steady-state rotations 

The equation of motion (^S)) admits solutions of the form which we will denote by ^i-g) (r) to distinguish 
them from the solutions 9{t) obtained in section 2] for F = 0. Fike in that section, we need only consider 
positive angular frequency solutions 0(g) (t, e, a), e > 0, describing direct counterclockwise rotations. Clock- 
wise rotations are then given by — 0(g)(r, — e, —a). For the purpose of perturbation theory we consider A, 
A and F to be small of the same order. In particular, F/A and A/A both remain finite as A — )■ 0. The 
perturbative computation is a simple extension of the analyses of sections l3 . 2 1 and \4A\ so we omit the details 
for brevity, limiting ourselves to quoting the results through second order. At zeroth order we have 

0(g)oir) = 9o{t) , (39) 

with 9q{t) given by (|27)) . so Qq is independent of F and relation ((28)) holds unchanged. In particular, for 
circular excitation there can be no contrarian motion, as in the case F = 0. At first order, 

^(s)iW = 0i(r)+F/(<,)i(r) , /(g)i(r)=sin(T + eo + «) , (40) 

with 9i{t) given in ([29]). At second order we get, 

OiM^) = 0(g)i(T) + \/^^h{T) + \^^Q2 + ^rV(g)2(r) + \T^Q(a)2 , (41a) 
with /2(t) and Q2 as defined in (|30|) and, 

n /^^tWn^ 1 1 - sin(e) sin(eo + 2a) 

e(,), = (^-tanOo) - 8 i + 3,,(,) ,,3(eo) 

/(g)2(T) = -^^(1 - sine) sin(T _ a) + ^(1 + sine) sin(T + a) (41b) 

5 A 1 
+ — — (1 - sine) sin(3T + 260 + a) + - sin(2T + 260 + 2a) . 

Notice that in this equation we have eliminated A in favour of Oq. 

A qualitative change with respect to the result of section |4] with F = is that the corrections f(g)n coming 
from gravity contain odd frequencies, so 9(^g){T) — r has period 2it. We notice also that the phase of 9(^g){T) 
acquires a dependence on F at second order. The perturbative solution to ([55]) described by ([55)) — ()^ 
exists only if (1^51) is satisfied, due to the definition (|?7)) of Bq- 
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5.2 Steady-state oscillations (F ~ A) 



We look for steady-state solutions of the form ([2^ . As in the previous section, we consider — 7r/2 < e < 7r/2 
thus including the case of circular excitation. For the purpose of the perturbative calculation we assume A, 
A, r < 1 to be small of the same order. From (l38t and (l22l) we find the zeroth-order solution 



^(5)o('r) = 0(g)o = "-TT - a , (42) 

with integer n. The value of 6(g)o results from consistency of the first-order equation, which then yields the 
first-order solution 

^(9)1 W = Q(ff)o + Ae(g)i + A/(g)i(r) , (43a) 
with /(g)i(T) given by (|33p and, from the requirement that f(g)2{T) be bounded, 

e(,),^-lf^Vcose)^ Bin(2e(,)o) 

4 ^ ^ cos(e(g)o + a) ^ ^ 

The second-order equation leads to 

0(g)2{T) = 9ig)iiT) + ^6(3)2 + ^/(<,)2(t) , (44a) 



with 



r 

2A 



0(g)2 = 77^^ — : — 7 ( ^sine- i(cos e)^ sin(eo - a) - ■^(cose)^ sin(3e(g)o 

cos['a(^g-^Q + a) \ L o o 

^^Q(9)i(cose)^cos(2e(g)o) 

/(g)2(7-) = (^6(3)1 - A/a) (1 + sine)cos(T - 8(3)0) + (o^g^i + X/Aj (1 - sine) cos(t + 8(3)0) (^^^^^ 

+ (j^^ ~ sin(T — a) — (1 + sin e) sin(r + a)) 

^(1 - sine) sin(T -f- 28(3)0 + «)-(! + sine) sin(T - 28(3)0 - a)) 

+ ^ ((1 - sine)' sin(2T + 28(3)0) - (1 + sine)' sin(2r - 28(3)0)) , 

where 8(3)2 results from requiring consistency of the third-order equation. 

In the case e = = a, corresponding to linear excitation in the same direction as gravity, ^(3)0(1") is an 
exact solution to (|38|) . the equilibrium solution. In that case the corrections 8(g)i_2, f{g}i,2 given in (j43|), 
(UD) vanish, as must happen to all higher-order corrections 8(3)^, f(^g)i with i > 2. For e ^ or a 7^ 
equations (|42]) - p4|) describe steady-state oscillations with dimensionless fundamental angular frequency 1 
(in physical units, angular frequency 7) about a central value 8(3) = Q(g)o + ^8(3)1 + A'/28(3)2 -I- 0{A^). 
For n even in (|42l) . the zeroth-order value 80 corresponds to the direction of gravity on the plane of the 
rotator. By numerically solving (j38l) we find that the solutions with n odd are unstable in this regime F '-^ A, 
as expected. ^ 

This perturbative solution is generally valid for A, A, F small, independently of whether condition (I28p 
is satisfied. If (p8)) holds, the steady state rotations of section [5TT] also exist and the system may reach one 
or the other steady state depending on its initial conditions. 



5.3 Steady-state oscillations (F <C A) 

The perturbative solution p2|) -(l44 )) given in the previous section was obtained under the assumption F = 
0(A). In particular, its limit as F ^> is not well defined. For F -C A perturbation theory must be 
accordingly modified, which we do in this section. We focus our discussion on the specific case A < 1, 
A = C'(A), F = 0{A'^). (Many other regimes are of course possible, which lead to results analogous to those 
obtained here.) With those assumptions for the parameters and (|22p for 6'(3)(r), we solve (|38p perturbatively. 
At zeroth order the solution is a constant, 0(3)o(t) — 8(3)0, to be determined. At order A, 

^(9)i(^) = 0(9)0 + A8(3)i + A/(3)i(r) , (45) 
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with /(g)i(T) given by and both 6(g)o,i stih undetermined. Consistency of the ©(A^) equation requires 
cos(e)2 sin(2eo) = 0, so that for e ^ ±t:/2 

n 

0(9)0 - , (46) 
for some integer n. With this, the second-order solution is 

^(9)2(t) = ^(3)i(t) + — (e(,)2 + /(g)2(r)) , 

f(g)2{T) ^-(^- 0(9)1^ (1 + sine) cos(t - 8(3)0) + + ^ia)!^ ~ s™*^) '^^^^'^ + ^(ff)o) ^4'') 
- ^(1 + sine)2 sin(2T - 26(3)0) + ^(1 - sine)^ sin(2r + 26(3)0) . 
The third-order equation requires 

^^^^^^ (co!e)^cos(2L(3)o) + ' 

for /(g)3 to be bounded. The expUcit expression for f{g)3 is rather lengthy, so we omit it for brevity. Finally, 
the ©(A^) equation leads to the consistency condition 

(=) 4 ^ 1 COs(6(3)o + «) 

^(^'^ = -^A^(cose)2 cos(26(3)o) ' ^^^^ 



Equations (|45p -(j49 l) yield a steady-state solution to ([38l) through second order in A < 1. Clearly, this 
perturbative solution breaks down as |e| approaches the upper bound Emax such that 6(g)i becomes ~ 1. 
Numerical study of these oscillating solutions suggests that for emax ^ |e| < 7''/2 they disappear from the 
spectrum of steady states. Therefore, as intuitively expected, the oscillations described in this section require 
an eccentric enough parametric excitation. 

In the case of linear excitation perpendicular to gravity {e — 0, a ^ ±7r/2) the central value 6(g) about 

which the pendulum oscillates was obtained in [TJ[20l[31] for A = 0. Our result 6(g) = 6(g)o + A6(g)i -I- 
A^/26(g)2 given by and agrees with those references at order A if we set A = in 

but further generalises it to the case A > 0, to second order in A, to any angle a, and to any large enough 
eccentricity (i.e., small enough e). The full solution (|45|) ~(|49 | . which gives the steady-state oscillations about 
the central value 6(g) through order A^, provides a similar generalisation to the oscillations described in fl8j 
for an inverted pendulum. From (|46|) we see that our results are applicable in the cases 6(g)„ — 2/c7r with 
integer k (normal pendulum) and 6(g)j, — {2k + l)n (inverted pendulum). As discussed in the next section, 
by numerically solving the solutions with 6(g)j, = (2fc -I- l)7r/2 are found to be unstable. 



5.4 Numerical results 

The steady-state rotations obtained perturbatively in section 15.11 arc compared with numerical solutions to 
([38ll in figures [3] (a) and (b). In figure [3] (a) we show 6{t) — t for linear motion of the suspension point 
(e = 0), at several angles a with the direction of gravity. The zeroth-order result (pQ)) is 60 — —0.025, 
independently of F and a. The first-order result (|40l) already accurately describes the numerical solution, 
with which the second-order result (|^T|) overlaps completely. In figure |3] (b) we show 9{t) ~ t for elliptic 
excitation with e = 7r/8 perpendicular to gravity, both for direct and contrarian motion. The zeroth-order 
result ((39|) is 60 — —0.018 for direct and -0.040 for contrarian motion. The first and second-order results 
are seen to be in tight agreement with the numerical solution. 

In figures [3] (c) and (d) we show steady-state oscillations as described in section 15.21 The parameters 
used in figure[2](c) do not violate the condition (1^ . but the initial angular velocity d9/dT{0) — 10^^ is too 
small for steady-state rotations (at e ~ tt/S). Numerically, for these parameters we find steady rotations for 
0.78 < d0/dT{O) < 1.44 for a = Btt/S, 0.75 < d0/dT{O) < 1.42 for a = 7/5, and 0.65 < d9/dT{0) < 1.55 
for a = tt/2. The rotator oscillates about a central position 6(g) whose zeroth-order value 60 — ~a is the 
direction of gravity, but receives small second-order corrections A^/26(g)2 = 10"'^, 6 x 10~^ and 10^^ for 
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a = 37r/8, 7/5 and 7r/2, respectively. If instead of e = 7r/8 we had chosen any other value — 7r/3 < e < 7r/4, 
the curves would vary by only a small amount. For |e| > 7r/4, however, steady rotations are obtained even 
for the small initial angular velocity used here. Notice that for e = ±7r/2 these oscillating solutions do not 
exist, so we expect that as e approaches those values the oscillating steady states should disappear from the 
spectrum of solutions. In figure [3] (d) we set the parameters to the same values as in figure [3] (c) except for 
a larger A, so that condition (1^51) is violated, and a larger initial angular velocity. Except for the buildup of 
additional turns during the transient, as indicated in the figure, the plots look almost the same as those of 
figure [3] (c). This is due to the fact that A enters the solutions of section [5^ onlv at order A^. 

The oscillating solutions for F ^ A obtained in section 15.31 are illustrated in figure HI We focus there 
on solutions with 6(g)o = (2n. + l)7r, describing the well-known phenomenon of the "inverted pendulum" 
(see the references cited in section [1] above) . In the case of linear parametric excitation parallel to gravity 
(e = = a), the exact solutions to ([55]) 0{t) = mi with n odd become stable for F <C A. In figure 2] (a) 
we show one such solution, including the transient and the relaxation to the steady state 0{t) = tt. When 
e = 7^ a there are no equilibria. In figure |4] (b) we show the oscillating steady state in that case, for 
various values of a. For e = and 6(g)o — the first-order result ()45p becomes constant, as shown in the 
figure, whereas the second-order result ()47|) accurately reproduces the numerical solution. The pendulum 
oscillates about 8(3) = 6(3)0 -I- ^^ig)! + AV2e(g)2 -h . . .~ tt -I- 1.2 x 10-^ tt -I- 2.3 x 10"^^ tt 2.9 x 10-^ 
TT-f 3.1 X 10^^, for a ~ 7r/8, 7r/4, 37r/8, 7r/2, respectively. We recall here that 9 is measured with respect to 
the trajectory of the suspension point, not that of gravity. For example, when a = 7r/2, figure |4](b) shows 
that the pendulum is oscillating about an almost horizontal position. Similar results are obtained in the case 
of elliptic excitation, as shown in figures |4] (c) and (d) with e = 7r/8 and e = 7r/4. In those figures we omitted 
the curves corresponding to a = 7r/2 for clarity, because they lie very close to those for a — 37r/8. As can 
be seen in the figure, the oscillation amplitude depends on both a and e. For e ^ 7r/4 higher orders beyond 
the second are needed to reproduce the numerical solutions and, as e approaches ±7r/2, perturbation theory 
becomes inapplicable. For e = ±7r/2 the oscillating steady states described in sections 15.31 do not exist. 

Finally, we notice that (|46p allows solutions oscillating about 6(g)o = (2n -I- 1)tt/2. Numerical study 

shows that, for the values of F, A, A considered in figure [H those solutions are unstable even for e = 0, 
and therefore they cannot describe steady states. We have also numerically solved for many sets of 
parameter values in the regime F <C A <C A, with e = and several angles a, and in all cases found the 
solutions to be unstable. That instability is of the same type as in the case F = 0, as shown in figure 
[H (d): the central value 8(g) about which ^(g)(T) oscillates is not constant, but its t dependence is slow 
compared to the period of oscillation. Thus, the oscillations with 6(g)o = (2n-|- l)7r/2 are well described by 
the perturbative solutions of section 15. 3[ but only up to an additive constant and for time intervals short 
compared with the time scales of variation of the central value 6(g). 

6 Parametric resonance 

Besides the steady-state solutions obtained in section [Sj when the value of F is in the vicinity of (^/2)^, t 
integer, (l38|) possesses resonant solutions. In this section we describe in detail the first parametric resonance 
mode, with F ~ 1/4. We assume A < 1, A = C(A) and F = 1/4 -I- 5 with 5 = 0(A), and expand in powers 
of A^/^ to find a solution of the form 



In this case it turns out to be more convenient not to separate explicitly from /fe(T) an additive constant term, 
as done in the expansions (fT2| . (|22|) . Substituting (|50|) into (|38|) . at zeroth order we find sin(8o + a) — 0. 
The zeroth-order solution is then. 




(50) 



k>l 



9o{t) = 9o = -a + mr , 
with integer n. At order A^/^ the equation of motion ([38|) yields 



(51) 




(52) 



dr2 
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Figure 3: Numerical solution (black solid line) to equation (p8| . compared with perturbative solutions at 
zeroth (grey dashed line), first (grey dotted line) and second (black dashed line) orders as given in section [5l 
Steady-state rotations (see sectionEl]): (a) e = 0, T = 0.1, A = 0.08, A = 0.001, 0{O) = 0.1, d0/dT{O) = 1.25; 
(b) e = ±7r/8, a = tt/2, T = 0.1, A 0.08, A = 0.001, 6(0) = 0.1, de/dT{0) = 1.25. Steady-state oscillations 
(see sectioning): (c) e = 7r/8, T = 0.1, A = 0.08, A = 0.001, 61(0) = 0.1, dd/dT{Q) = lO^^ (d) e = 7r/8, 
r = 0.1, A = 0.08, A = 0.06, 61(0) = 0.1, d0/dT{O) = 10. 



Requiring /i to be bounded leads to the condition that n in (|5ip must be even, so the pendulum oscillates 
at resonance about its position of lowest potential energy. The first-order solution is then, 

0i(r)==eo + Ai/Vi(r) , /i(t)-picos(t/2-(pi) , (53) 

with the constants pi, ipi undetermined at this order. Similarly, the equation of motion at order A is 

^Jl + lf^{T) = -i(l-sine)sin(T-a) + i(l + sine)sin(r + a) , (54) 
leading to the second-order solution 

92{t) = 0i{t) + Ah{T) , 

2 2 (55) 

/2(t) = P2 cos(t/2 - </'2) + ^(1 - sine) sin(T - a) - -(1 + sine) sin(r + a) , 

with the constants p2, '■P2 undetermined at this order. At ©(A'^/^) we obtain the equation for /3(t) 



d7^ + 4^^(") = -Ad7~A^^(") + 24^^^") 2 



h{r) = --^^ -Mt) + -MtY - -(1 - sine)cos(T - a)/i(r) 

^ ' (56) 
- -(l + sine)cos(r + a)/i(T) . 

The right-hand side of this equation, with /i(t) given by ()53p . contains terms proportional to cos(r/2) and 
sin(T/2) leading to an unbounded solution fsir). Requiring those terms to vanish determines the constants 
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Figure 4: Inverted-pendulum solutions to equation ([38]) . Numerical solution (black solid line), compared 
with perturbative solutions at zeroth (grey dashed line), first (grey dotted line) and second (black dashed 
line) orders as given in section [5T3l The parameters are set to F = 10"^, A = 0.08, A = 0.0375, 9{0) = tt, 
d0/dT{O) = 10-^ and (a) e = 0, a 0; (b) e = 0; (c) e = tt/S; (d) e = 7r/4. 



pi, ifii in (|53l) to be 



pi = 4:\/2S/A + y (cosa)2 + (sine)2(sina)2 - (A/A)2 , 

- cos a + y (cos q;)2 + (sine)2(sina)2 ~ (A/A)^ 

(fix = arctan | — I + ti^^tt , 

A/ A — sine sin a 



(57) 



where n^-^ is an integer parameter on which the solution depends, similarly to the parameter n in (j51|) . We 
remark, however, that from (j53l) it is apparent that the solutions depend only on whether Uip^ is even or 
odd. With (|57|) . equation leads to a third-order solution of the form. 



^3(r)=02(T)+A3/2/3(r) 



1 3 ...... 1 



/3(r) = A3 sin(T/2) -1- B3 cos(t/2) - —p cos(3r/2 - 3<^i) + -p{l - sine) cos(3t/2 - a - (fi) (53) 

-I- -p{l + sine) cos(3t/2 + q — ipi) . 
8 

The coefficients A3, B3 in /3(t) are determined by the requirement that the perturbative equation of order 
/\5/2 jjj^ygi; Jiave a bounded solution. The resulting expressions are too lengthy to transcribe here, so we 
omit them (see, however. Appendix]^. Finally, consistency of the equation at order A^, which determines 
/4(t), requires 

P2 = , (59) 

which fixes the second-order solution 92(t) in (|55p . Summarising, the steady-state solution at order A° is 
given by dsn; at order A^/^ by ^ with ^] at order A^ by ^ with ([591); and at order A^/^ by 
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We remark that /i(t) oscillates with angular frequency 1/2, and f2{T) with angular frequency 1. /i,2(''") 
determine the equations for the higher-order corrections fkir), k > 2. The structure of those perturbative 
equations is such that /fe(r) contains only terms with half-integer (respectively integer) angular frequency 
if k is odd (respectively even). As a consequence, the following relations, which are easily checked for the 
low-order solutions given above, actually persist to all orders of perturbation theory 

6{t, e, a, n^^ + 1) = -d{T, — e, —a, n^J + 2mT — 6{t + 2tt, e, a, n^J . (60) 

From the second equality we immediately see that in the case of linear motion of the suspension point parallel 
to gravity, e = = a, all terms with integer angular frequencies must vanish. In that case only the terms 
with odd k contribute to (ISUl) . In general, for arbitrary values of e and a, equations describe two 

different solutions according to whether is even or odd. Those two solutions are related, modulo 27r, by 
the symmetries of as shown by (pD|) . Notice that, whereas the invariance of under 9{t) 0{t + 2-k) 
did not play any role in section [5] because there only integer angular frequencies appear, that symmetry acts 
non-trivially on the resonant solutions described in this section. 

We see from ([Ff|) that there are two necessary conditions for these perturbative solutions to exist. The 
first one places an upper bound on damping for a given excitation amplitude, 

A < AV(cosa)2 + (sine)2(sina)2 . (61) 

Thus, for linear excitation (e = 0) there is no parametric resonance if a = ±7r/2. For circular excitation the 
right-hand side of (j61l) becomes independent of a, as expected. The second condition, 

5 > -iAy^(cosa)2 + (sine)2(sina)2 - (A/A)2 , (62) 

puts a (negative) lower bound on S — T — I / 4. Notice that there is, apparently, no restriction on how large S 
can be as long as it is positive. This is, in fact, an artifact of perturbation theory: a non-perturbative stability 
analysis leads to well defined limits on the magnitude of S for resonance to be possible [21]. Yet, perturbation 
theory does restrict how negative S can be through ()62p . Numerical study confirms this asymmetric situation: 
given < Sq < A such that — violates (|62p. resonant solutions can be found numerically for 5 = Jq, but 
not for S = —Sq. 



6.1 Steady-state rotations with angular velocity 1/2 

When r = 1/4 -f ^, with 6 = 0(A), A < A and A <C 1, equation (|38|) can be seen as a small perturbation 
about the unperturbed equation 

— + uj^ sm{e{T)) = , a; = 1/2. (63) 

In that case there must be solutions to (pS)) close to those of with corrections of 0{A). For oj^ ~ 1/4, 
possesses solutions of the form 6{t) = t/2 + Qq -|-sin(r/2 + Qq) + . . ., for some constant Qq, that can be 
computed perturbatively as an expansion in powers of lo'^ . How the remaining terms in psp are added to (I63p 
as perturbations depends on their relative size as compared to a;2 = 1/4. For A ^ 10^2-io~i^ A ^ 10"'^- 
10~2j (5 ^ A, an appropriate expansion is obtained by setting A = 0{uj'^), 6 = ©(w^) and A = 0{uj^). Thus, 
we expand psp perturbatively, in a completely analogous way as in the previous sections, with 

oo 

0{r) = -T + eo + Y^i^^riQn + Mr)) . (64) 

^ n=l 

At zeroth order we obtain, 

6'o(t") = ^T-f Oo , Qo ^ -^arcsin ^ ^^^ + sine) j ~ ^'^ + ' (6^) 

with n integer. We see that these steady states exist only if 

A<A(H-sine). (66) 
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At first order we get, 



0i(t) =0o(r) + -(ei+/i(T)) , -ei = -2(5tan(2eo + a) , -/i(r) 



sin(T/2 + eo + a) . 



(67) 



We remark that this first-order solution aheady depends on all of the parameters in p8|) . and therefore on 
all of the physical quantities entering it: elliptic parametric excitation and gravity, their relative angle, and 
damping. Notice also that at this order (l67l) is a solution to ((63|) . the additional terms in (l38l) entering (|67| 
only through 0o,i- At second order the perturbative solution is 



Because the expansion parameter = 1/4 is not very small, the third-order correction is necessary in order 
to obtain a perturbative solution as accurate as those given in previous sections. We take into account the 
third-order terms in numerical computations in the following section, but we omit them here for brevity. 

6.2 Numerical results 

The previous results on parametric resonance are compared to numerical solutions to (|38)) in figure [5j The 
case of linear excitation parallel to gravity, e = = a, is shown in figure[5](a), where we plot the order A^/^ 
result (|53|) . with ([57l) . and the order A'^/^ result given in (f58| and Appendix [Al The parameters in that 
figure were chosen so S is very close to the bound (p^ which in this case reduces to 5 = —0.047 > —0.05. 
This is the region where perturbation theory begins to break down, thus appropriate to test the 0{A'^^^) 
correction. As seen in the figure, the numerical solution is accurately reproduced by our perturbative result. 
For slightly higher values, S > —0.04, the 0(A^/^) result would be enough to achieve the same accuracy. 

In figure [5] (b) we illustrate the case of linear excitation not parallel to gravity, e = ^ a, for which the 
order A^ correction (I55p . with ([59|) . is non- vanishing. The agreement of the perturbative results with the 
numerical one is very accurate. The decreasing amplitude of oscillation with increasing a is apparent from 
the figure. As discussed above, (j6ip implies that there is no resonance when gravity is perpendicular to the 
linear trajectory of the suspension point. This is confirmed in the figure, where the numerical solution with 
a = 7r/2 is seen to be an oscillating solution of the type described in section [521 with twice the frequency 
of the resonant solutions and a much smaller amplitude. 

The case e 7^ is shown in figure [5] (c), where we set e = 7r/2 corresponding to counterclockwise 
circular motion of the suspension point. Unlike figure [5] (b), the oscillation amplitude does not show a 
discernible dependence with a, and for a = 7r/2 a resonant solution is obtained. Aside from those quantitative 
differences, figures (b) and (c) are qualitatively similar. The behaviour for other values of e is intermediate 
between those two. 

In figures [S] (d), (e) and (f) we show steady-state rotations with angular frequency 1/2, as described in 
section [Q] It is clear from the figures that the dependence on e and a is rather weak. The lower bound dHU) 
is violated in figures (d) and (f), so that a resonant steady state is not possible in those cases, but it is not 
violated in (e), which shows that steady states with angular velocity 1/2 and resonant ones can coexist. If we 
had chosen d9/dT{0) slightly different from 1 in figure (e), the solution we would have obtained numerically 
would have been a resonant one. In fact, finding these steady states numerically seems to require fine tuning 
of the initial conditions, that being the reason why we chose the same ones in the three figures. 

7 Final remarks 

In the foregoing sections we discussed several steady-state solutions to the equation of motion (j4|) with 
parametric excitation given by ([3]). We restricted our treatment to the parametrically excited rotator. 




e2(r)-0i(r) + -(e2 + /2(r)) , 



(68) 



- sin(r + 260 + 2a) - 2A(1 + sine) sin(r/2 - Oq) + -A(l - sine) sin(3/2T + Gq) 
8 9 

+ 45sin(T/2 -f 60 + a) -I- ^61 cos(t/2 + 60 + 0;) . 
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Figure 5: Solutions to equation (j38]) in the parametric- resonance region. Numerical solution (black solid 
line), compared with perturbative solutions at first (grey dashed line), second (grey dotted line) and third 
(black dashed line) orders as given in section [51 The parameters are set to (a) e = 0, a = 0, n^j — 0, 
r = 0.203, A = 0.1, A = 5 X lO^^, 9(0) ^ 0, de/dT{0) = -10-*; (b) e = 0, n^, = 0, T = 0.253, A = 0.01, 
A = 0.005, 9(0) = 0, de/dT{0) = 0.01; (c) e = 7r/2, n^, = 1, T = 0.253, A = 0.01, A = 0.005, 6(0) = 0, 
de/dT{0) = 0.1; (d) e = 0, a = O^F = 0.243, A = 0.01, A = 5 x lO^^^, ^(O) = 0, de/dT{0) = 1; (e) e = 0, 
a = 7r/8, r = 0.253, A 0.01, A = 5 x 10-^ ^(O) = 0, d9/dT{0) = 1; (f) e = 7r/8, a = 0, T = 0.244, 
A = 0.01, A = 5 X 10-^ 61(0) = 0, d6i/dr(0) = 1. 
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(A < 1 in dH)), as opposed to the orbitally excited rotator (A » 1), and to low damping, 10~^ ^ ^ ^ 1- 
Similarly, for the pendulum we considered also F < 1. As mentioned in section [U many other regimes 
are possible. Some of those may be treated as simple variants of the perturbative methods used here, as 
occurs, for example, with the two versions of the condition F < 1 given in sections 15.21 (F — 0{A)) and 15.31 
(F = 0(A3)), which lead to slightly different perturbative expansions, with qualitatively different results. 
Yet, other parameter regions not considered here may require other mathematical tools. But even within 
the domain of parameter space considered here, and even within the class of simple steady states with 9{t) 
either monotonic or periodic, our results are certainly not exhaustive. 

In section |3] we discussed the linearly excited parametric rotator. We gave the necessary condition for 
steady-state rotations and showed that 9{t) — r is a periodic function with period tt that oscillates about 
a value O determined by the parametric amplitude and the damping coefhcient. There is also in that case 
a steady state in which the rotator oscillates about the transverse position. The latter solution seems to 
be unstable in the presence of gravity or eccentricity. In section 2] we extended those results to the case 
of elliptic motion of the suspension point. For steady-state rotations we generalised the results for linear 
excitation to any eccentricity, including circular excitation, for direct motion, and established the existence 
of contrarian steady state rotations for any elliptical excitation except the circular one. We treated also 
steady-state oscillations with the same angular frequency as the parametric excitation. Those oscillating 
steady-state solutions are characteristic of strictly elliptic excitation (0 < |e| < n/2), since they do not 
exist for linear or circular motion of the suspension point. In the linear case the oscillating solutions reduce 
to equilibrium ones. The case of circular parametric excitation possesses some special characteristics that 
cannot be treated with the perturbative method discussed here, so we will discuss it separately [55] . 

In section [5] we considered the parametric pendulum. For steady-state rotations the perturbative theory 
with F ~ A is only a slight modification of that with F = of section |31 with corrections due to gravity 
inducing odd angular frequencies in the spectrum of oscillations about uniform rotation. We obtained 
explicit expressions for those steady states, as functions of the eccentricity and angle with gravity. As in the 
case without gravity, contrarian-motion solutions exist for all e and a, except for circular excitation. The 
steady-state oscillatory solutions with F ~ A are very different from the case F = 0, because &o in that case 
is completely determined by gravity. By contrast, in the case F ^ A the oscillating solutions are a small 
perturbation about the case F = and, in particular, Qq is independent of gravity. This is the essential 
feature of the inverted-pendulum phenomenon: because F ^ A, the dynamics are mainly determined by 
the torques of the inertial forces due to the motion of the suspension point, and not by those of gravity. 
We compute the angle 6 about which these steady states oscillate as a function of the eccentricity of the 
suspension point's trajectory, as measured by the parameter e, and its angle a with gravity. That result 
provides an analytical basis to the observation that, as a not-yet-inverted pendulum is slowly rotated from 
the downward to the upright position, it follows the motion of the parametric excitation with an angular 
lag [18]. More generally, as discussed in quantitative detail in that section and is intuitively expected, the 
dependence of the solutions to Q on the angle a between gravity and the major axis of the elliptic excitation 
is stronger for linear parametric excitation, and grows weaker as |e| is increased. 

In section[6]we considered parametric resonance, for elliptic motion of the suspension point with arbitrary 
eccentricity and at an arbitrary angle a with gravity. We found the explicit form of the oscillatory steady 
states in the resonance region, and their dependence with eccentricity and a. Unlike the usually studied 
case of linear excitation parallel to gravity, when the excitation is not linear or not parallel to gravity, 
integer angular frequencies enter the spectrum of oscillations, besides the half-integer frequencies familiar 
from the case e — = a. We obtained necessary conditions for parametric resonance in terms of e and a. 
In particular, for linear excitation the amplitude for parametric resonance vanishes at a = 7r/2, but not for 
other eccentricities. We found also a steady-state rotation mode with angular frequency equal to half the 
parametric excitation frequency, in the parameter region associated with parametric resonance. Unlike the 
resonant solutions, however, the dependence of the steady-state rotations described in section 16.11 on the 
parametric excitation amplitude is analytic. Those rotation modes do not seem to have been discussed in 
the previous literature. 

Finally, we want to stress that the steady states discussed above are only a handful out of a much more 
numerous spectrum. For linear and highly eccentric elliptic excitations there are other modes, even at small 
parametric excitation, that are of a rather different nature than those considered here and will be discussed 
separately. For circular and slightly eccentric elliptic excitations, as mentioned above, an extension of the 
methods presented here is necessary to uncover some of the simplest basic motions [36] . On the other hand. 
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for the steady states considered in this paper, the regions of stabihty and basins of attraction were probed 
rather crudely by means of numerical solutions to the equation of motion. More detailed and technically 
accurate analyses are required in those respects, which lie beyond the scope of the present paper and which 
we defer to future work. 
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A Parametric resonance with e = = a 

In section [6] a solution 9(t) at parametric resonance is given through 0{A). For e = = a, however, the 
terms of 0{A) vanish, as explained in that section, so the correction at 0(A'^/^) becomes important. In this 
appendix we explicitly provide the coefficients A3, B3 in (|58p in that case. 
When e = = a, consistency of the ©(A^/^) equation for /5(t) requires, 

Na = bi cos fi + fli sin</?i + 63 cos(3(^i) + sin(3(pi) , 
bi = -16/ipi(256 - lUpl - 256dpl + 3pt) , 

ai = pi(-2i2 _ 2i3^ _ 2iip2 _ 2^ dpi + 2^^d^pl + 192pt - 352dpf + ipl) , 
63 = 768V? , 03 = 48p? (16 + 32d + pi) , 
Nb = b'l cos ipi -\- a'l sin ipi + 63 cos(3(^i ) + ag sin(3(,ci ) , 
b[ = pi(2i2 - - 2"p? + 2^dp? + 2i3rf2p2 _ ^Q2pj - 352dpj + 3pf) , 
a'l = 16Vi(256 + lUpl - 256dpl + Spf) , 
b'^ = 48p? (-16 + 32d + pI) , 4 = -768V? , 
D = 64(-2^ + 2i°d2 + 2^h^ - 2'' dpi + ip\ + Z2pl cos(2(^i)) , 

with pi, Lpi as defined in (|57p . and with h — A/A, d = (5/A. 
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